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PREFACE 



This report was written to supplement material by 
S. M. Pollock in Selected Methods and Models in Military 
Operations Research as well as material in OEG Report No. 56 by 
B. 0. Koopman. Much of the report is based on a course given 
by S. M. Pollock at the Naval Postgraduate School in 1969. 

The report's relation to OEG Report No. 56 will be evident on 
reference to that work. 

Some other sources of material to which this report 
relates are listed in the Bibliography. 
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I. Detection Theory and Detection Models 

In signal detection theory, the decision making portion 
of a detection system is called the receiver and a detection 
experiment is the observation by a receiver of its input during a 
time interval. The input which is related to a target is called 
signal, and the input which is not related to the target is called 
noise. In general, the observation is assumed to be of a known 
region which in some cases is called a resolution cell. 

Signal detection theory is a basis for detection modeling. 

The detection models which are discussed here rely heavily upon 
it. -When a detection experiment has been performed either the 

event = {At least one target was present in the region which 

was observed during the time of the observation.} or its comple- 
ment Hq will have occurred. In general, the models will specify 
that either the event D^ = {The receiver decides at least one 

target was present in the region which was observed during the 

time of the observation.) or its complement Dq will have occurred. 
In terms of signal and noise, the events Dq and D^ can be 
expressed as follows: Dq = {The receiver decides its input during 

the time of the observation was noise.} and D^ = {The receiver 
decides its input during the time of the observation was signal 
and noise.} The theory which is the basis for detection models 
which have the above properties is called binary detection theory. 

Four events which are important in binary detection theory 
are indicated in the Venn diagram of Figure 1. 




Figure 1. Four events important to binary detection theory. 

The Venn diagram emphasizes a decision problem which is 
associated with a receiver. The problem is this: Under what con- 
ditions should the event occur, that is, under what conditions 

should the receiver decide that a target was present during the 
time of an observation? Detection theory may be able to provide a 
solution to this problem. If this is the case, the noise will 
generally be described as a random process and the signal will be 
described as either a deterministic or a random process. 

The following notation and terminology will be used: 

Pf = P(D 1 | H q ) will be called the false alarm probability, 
p^ = P(DjJh^) will be called the detection probability. P = P(H^) 
will be called the prior probability. P is the probability that 
a target will be in the region to be observed during the time of 
the observation. 

The input to a receiver is assumed to be a quantity whose 
square is proportional to power. The input at some time t^ will 
be symbolized by y(t^). The noise at t^ will be symbolized by 
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n(t^) and the signal by s(t^). The noise and signal processes 
will be represented by {n(t i ),t i et} and {s(t.),t.et} where t 
represents the time interval during which the receiver 
observes the input. Often s(t^) and n(t^) can be assumed to 
be additive in which case H Q = {y(t^) = n(t^),t^et} and 

H 1 = {y(t i } = + s(t ± ) ,t i et). 

Because of the finite quantity of information present at 
the input of a receiver, y needs to be measured at only a finite 
number of points in time in order to be adequately determined over 
an observation interval. For this reason, the noise and signal 
processes can be represented by (n(t^), i = l,...,m) and 
{s(t^), i = l,...,m) where t^,...,t m are in the observation 
interval of length t. The noise process is then defined by a set 
of m random variables. To specify the noise process, one needs 
to specify only the joint distribution of the m random variables. 
If the signal process is deterministic, it is simply a finite set 
of values which is known before the detection experiment is per- 
formed. 
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II. Decision Criteria 

To simplify the discussion of decision criteria and decision 
rules, the stochastic process which represents a receiver's input 
will be assumed initially to be a single random variable Y. It 
can be called the decision random variable. The input process in 
this case is determined by the two conditional distribution func- 
tions F y (y | Hq) and F y (y|H 1 ). 

The condition that a receiver's input is required to satisfy 
in order that will occur can be specified in terms of a deci- 

sion rule. For the assumed case, a decision rule is a rule which 
determines for every observable value of Y the decision that the 
receiver is to make. The decision rule can be considered to be a 
function <p (y) which relates the observable values y to the 
following two statements: 

dp.* Decide the input was noise. 

d^: Decide the input was signal and noise. 

Choosing a decision rule <p (y) defines a set £2 such that 
D ± = {Y e £2}. 

The problem which was considered above can now be restated 
in the following way: What criterion should be adopted in order 

to determine a decision rule? A desirable characteristic for a 
criterion is suggested by the following argument: Consider the 

odds in favor of given y is observed. That is, consider 

P(H 1 |y)/P(H 0 |y) . 
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One might expect that the values of y in £2 should make this 
ratio relatively large. But making this ratio large is equiva- 
lent to making the likelihood ratio L(y) large. This suggests 
that ft might be defined as follows: ft = {y: L(y) ^ K) where 
K is some positive constant which has yet to be specified. 

Four specific decision criteria are defined next in terms 
of ft. For each criterion ft has the above form with only the 
procedure for determining K being different. The decision cri- 
teria are: 

1. The Neyman-Pearson Criterion: Choose ft so that p^ is 

a maximum subject to the constraint that p^ ^ a where a is 

a specified value. For a continuous decision random variable, 
the constant K is chosen so that p^ = a. 

2. The Bayes Criterion: Choose ft so that the expected cost 

of a receiver's decision is a minimum. For a continuous decision 

random variable, K = ( c ^o -C 00 //C 01 -C ll^ if c 1Q > c Q0 and 

C Q1 > C.^ where C^j is the cost of the event fi H ^ . 

3. The Ideal Observer Criterion: Choose ft so that the 

probability that the receiver makes an incorrect decision is a 
minimum. For a continuous decision random variable, K = (1-P)/P. 

4. The Minimax Criterion: Choose ft so that the maximum 

expected cost of a receiver's decision is a minimum. For a con- 
tinuous decision random variable, K = ( c io -C 00 //< '01~ C ll^‘*' -P * ^ 

if C 1Q > Cqq and C Q1 > C^. Here P* is the value of P which 
would make the expected cost of a receiver's decision a maximum 
if a Bayes decision rule cp B were used, and P is unknown. 
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A more general discussion of the above criteria would involve 
the notion of a randomizing rule. 

If a model is adopted which specifies the conditional dis- 
tributions and a decision rule, then the values of p f and p^ 
can be determined in principle . The pair of values (p^,p^) is 
called a receiver operating point. If the decision rule results 
from one of the four likelihood ratio criterion listed above, then 
it will involve the parameter K through the relation 

ft = {y: L(y) £K}. And, for a given value of K, since ft determines 
(Pf'P^)/ a single operating point results. By varying K, a set 
of operating points can be generated which is called an ROC curve. 
Different ROC curves can be produced by changing one or both of 
the conditional distributions. Note, a change in a conditional 
distribution function implies a change in the signal or the noise. 

Decision rules which result from applying the four likelihood 
ratio criteria in a model in which the input process is determined 
by a set of m random variables can be expressed in terms of a set ft 
as follows : ft = { (y^ , . . . , y m ) : L (y^ , . . . ,y m ) ^ K} where K is 

specified in the same way as it is in the corresponding case in 
which the input process is determined by a single random variable. 
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III. Two Detection Models 

A particular detection model will now be examined. In 
the model, the receiver's input process is defined by: 

Hq : {y (t^) = n (t^) , i = l,...,m} and 

H l : {y (tj.) = n + s ( fc i) / 1 = 1/ • • • fin) 

where the values of y are measured every At units of time. 

The noise process {n(t^), i = l,...,m} is assumed to consist 
of a set of m independent normal random variables each with 
mean zero and variance a 2 . The signal process {s(t^), i = l,...,m} 
is assumed to be deterministic. Thus, the input process consists 
of m independent normal random variables Y^,...,Y each with 
variance a 2 and each with mean zero when the target is not 
present and with mean s^ = s(t^) when the target is present. Note, 

such a model might be used to obtain an optimistic estimate of a 

detection systems performance, since all the information about 
the signal is assumed to be known. 

The result of the application of a likelihood ratio decision 

rule in the model can be expressed in terms of a random variable 

Z. This random variable is called the cross correlation statistic 

and it is defined by Z = £ s i Y i* However, it is more convenient 

to express the results in terms of a random variable V which is 

defined by V = Z/a . In terms of this random variable, 

z 



and 



p f = 1 - $>(v*) 
P d = 1 - $(v*-d%) 
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where $ symbolizes the standard normal distribution function and 
v* = (1/0 ) (0 2 In K+ (1/2) £ s 2 ) , d ^ = I s?/o . 

Z X X 2 

The parameter d is called the detection index and, since 

°l = l s i° 2 ' d = l s [/o 2 . 

Recall, the input stochastic process is assumed to repre- 
sent a quantity whose square is proportional to the input power. 
Therefore, the average power input to the receiver over the time 
interval t during which the receiver observes the input from a 
region can be approximated by £ y?/m. The average signal power 
can be approximated by S = £ s?/m and the average noise power 
by I n?/m. And the expected average noise power can be approxi- 
mated by N = l 0 2 /m = 0 2 . In these terms, d = m(S/N) . 

The receiver's bandwidth will be represented by BW. 

Assuming m corresponds to the value determined by the sampling 
theorem, that is, m = t/At = 2t (BW) , the detection index can be 
written as d = 2t(BW)(S/N) or, defining Nq = N/Bw, as d = 2t(S/Ng). 

A model in which the signal process is not deterministic 
will be considered next. The model is defined by 

V {y i = n i' i = 1, . . . ,m} 

H 1 *. {y i = n i + s i , i = l,...,m}. 

Here the noise process is a set of m independent normal random 

variables with mean zero and variance 0 2 as before. However, 

now the signal process is also a set of m independent normal 

random variables with mean zero but with variance o 2 . The result 

s 

of applying a likelihood ratio decision rule in this model can be 
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expressed in terras of a random variable X which is defined by 

X = £ Y?. When a target is not present, the statistic X/a 2 has 

a chi-square distribution with m degrees of freedom and when a 

target is present, the statistic X/(o 2 +o 2 ) has a chi-square 

s 

distribution with m degrees of freedom. Hence, p^ = P(x^x*/a 2 ) 

and p, = P [x 2 ^x*/ (c 2 +o 2 ) ] where x 2 is a chi-square random 
a m s J m 

variable with m degrees of freedom and x* is a number which 

is determined by the decision rule. 

When the target is not present, since X/o 2 has a chi-square 

distribution, the variance of X is 2mo 4 and the mean is mo 2 . 

When the target is present, since X/(o 2 +o 2 ) has a chi-square 

s 

distribution, the variance of X is 2m(o 2 +o 2 ) 2 and the mean is 

s 

m(o 2 +o 2 ). By the central limit theorem, as the number of degrees 
of freedom m of a chi-square random variable becomes large it 
can be approximated by a normal random variable with the same mean 
and the same variance. Hence, for sufficiently large m, p^ 
and p^ can be approximated as follows: 

p f =; 1 - $ [ (x*-ma 2 ) / (2ma 4 ) 

and 

P d = 1 - ${ [ (x*-mo 2 ) - mo|] / (2mo 4 ) ^'} 

where (o 2 +o 2 ) 2 has been approximated by a 4 . This approximation 
can also be written as: 

p f = 1 - $ (v*) 

P d = 1 - $ (v*-d' 5 ) 

with v* = (x*-mo 2 )/ (2mo 4 ) ^ and d = (mo 2 ) 2 /2mo 4 . 

b 
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Note, over the observation interval t, £ s?/m 

the average signal power so S=£a 2 /m=o 2 is 

s s 

the expected average signal power and, as before 
approximately the expected average noise power, 
model the detection index can be written as 

d = j (S/N ) 2 = t (BW) (S/N) 2 . 



is approximately 
approximately 
N = a 2 is 
Hence, for this 
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IV. Detection Model Applications 

When a likelihood ratio decision rule was used with either 
of the two models discussed above, the following result was obtained 

p f = 1 - $ (v*) 

p^ = 1 - $ (v*-d 2 ) 

In the case of a known signal (deterministic process) , a cross 
correlation receiver is required, and d = 2t (BW) (S/N) . This will 
be called Case I. In the case of a Gaussian signal, a square law 
receiver is required and, for S/N << 1 and t(BW) large, 
d = t(BW) (S/N) 2 . This will be called Case II. In both Case I 
and Case II, the detection index d is a function of (S/N). 

The two models can be related to radar and sonar systems, for 
example, by using the radar equation or sonar equation to deter- 
mine (S/N) . 

In some radar and sonar detection models, it is assumed 
that the decision rule used results in a required p^ and in 
addition that a minimum value of p d is specified in the sense 
that a target is said to be detectable only when p^ is greater 
than or equal to the minimum value. This minimum value of p^ 
along with the required p^ define what could be called a mini- 
mum acceptable signal-to-noise ratio (S/N) m . In sonar models, this 
minimum acceptable signal-to-noise ratio determines the detection 
threshold DT by the relation DT = 10 log(S/N) m . If the mini- 
mum acceptable value of p^ is .5, then DT is usually called 
the recognition differential RD. 
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A passive sonar detection model will now be discussed. It 
is an example of a model which uses the concept of a minimum value 
of p^. It also illustrates a method of dealing with non-station- 
ary noise and signal processes. In the model, the event detection 
is the event {Xg E ^0} where X gE represents signal excess, a 
random variable with expected value SE. The source level, noise 
level, directivity index, and the detection threshold or recogni- 
tion differential are also random variables with expected values 
SL, NL, DI, and DT respectively. The passive sonar equation 
can be written in terms of these random variables or in terms of 
their expected values. In the latter terms, SE = FOM - TL, where 
FOM = SL - (NL-DI) - DT is the figure of merit. To determine the 
probability of a detection PCX.^O), the distribution of the 

random variable X must be specified. 

bb 

This model could be considered to be equivalent to one in 
which the event {X SE ^0} is the event that p^ is greater than 
a minimum acceptable value. Although this definition may appear 
more reasonable than the one which defines {X.^O} as the event 
detection, it would likely be no more useful for developing tac- 
tics, designing exercises, or testing detection devices. 

The noise and signal processes should be effectively sta- 
tionary during a "look" at a resolution cell, that is, over an 
integration period if a Case I or Case II model is used to relate 
Pf, p^, and the signal-to-noise ratio. If either of these models 
is used, the signal-to-noise ratio for an integration period will 
be the value of a random variable, and, therefore, p^ will be 
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conditioned on the random variable signal-to-noise ratio. The 



receiver is usually assumed to be readjusted so that 
relatively constant from "look" to "look" at a given 
cell in applications of the model. 



p^ remains 
resolution 
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V. Search Theory and Search Models 

Search theory deals with problems such as that of determining 
the optimal allocation of search effort in a region or that of 
determining the probability of detecting a target within a given 
length of time for a particular search plan. In treating these 
problems, false alarms are often ignored. In effect, the time 
spent on false alarms is assumed to be small relative to the time 
spent on the search, however, the cost associated with a false 
alarm is not assumed to be negligible. 

In general, it is assumed here that during a search a 
target will move through resolution cells surveyed by the detec- 
tion system of a searcher. In the model of such a search, the 
probability that the input from a resolution cell containing the 
target will cause the receiver to decide the target is present 
is p^. In some detection systems, the receiver can delay the 
decision regarding the presence or absence of a target in a reso- 
lution cell and can recall the input of adjacent resolution cells. 
Such detection systems are difficult to model. To illustrate 
this, consider the following hypothetical system consisting of an 
active sonar system plus an operator. Assume that p d as a func- 
tion of the signal-to-noise ratio has been determined for the 
operator in a laboratory experiment by forcing the operator to 
respond after a single "look" at a resolution cell. Operationally, 
the probability that the operator detects a target will not be 
given by the laboratory determined value of p^ if the operator 
delays the detection decision for several "looks." Various models 
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have been proposed to deal with this situation. In one, a "three- 
out-of-five" detection criterion has been adopted. The criterion 
states that an operator will declare a target if out of five con- 
secutive inputs from adjacent resolution cells at least three 
are such that they would have caused the operator to declare a 
target present in a single "look" (forced choice) experiment. An 
input which satisfies this condition is referred to as a success. 
With this criterion, the probability that a target will be detected 
with m consecutive "looks" at the target is found as follows: 
Determine the 2 m sequences of successes and failures of length 
m which could result with m consecutive "looks." The proba- 
bility is equal to the sum of the probabilities of occurrence 
of those sequences for which the three-out-of-five criterion is 
satisfied. The probability of a particular sequence will depend 
on the probability of a success and the probability of a failure 
for each of the m "looks" in the sequence. Unfortunately, in 
most situations, the number of sequences to be considered is too 
large to be computationally tractable even if the probabilities 
of success and failure for each "look" in the sequence can be 
determined . 
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VI. The Probability of Detection During a Search 

The problem of determining the probability that a target 
will be detected in a search at or before n "looks" have occurred 
is generally basic to the solution of search problems. Let N 
represent the number of the "look" at which detection first occurs, 
then this probability can be expressed as 

n 

P(N£n) = l P (N=i) . 
i=l 



It can also be expressed as 



n 

P(Ns£n) = 1 - TT (1-g, ) 
i=l 

where g^ = P (N=i | N^i-l) is the probability that the event detec- 

th 

tion occurs on the i— look conditioned on the event detection 
has not occurred earlier. The second expression for P(N^n) is 
generally of greater interest than the first one, since g^ can 
be directly related to operational parameters such as target range. 

j_ 'U. 

Note, if the resolution cell being examined on the i— look is 
empty, then g^ will be zero. 

A continuous analog to P(N^n) can be developed as follows: 

Consider a model in which the time for a look is At, and if 

th 

detection occurs on the i— "look," the time of detection is 

4- Vi 

iAt. The probability of detection at or before the n— "look" then 
can be expressed in terms of time by noting that P(N^n) = P(T^nAt) 
where T represents the time of detection. Now define the detec- 
tion rate to be y (nAt) = (1/At) P [T=nAt | T^ (n-1) At] . The probability 
P(T^nAt) can be approximated in terms of a related quantity y(t) by 
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P (T^t) = 1 - exp{- 



t 



Y (T)dT) 

0 

where Y(t) is related to Y (nAt) by the requirement that 
P(T£t) equals P (N^n) for t = nAt. In some models where it 
would be useful to use P(T^t), the continuous approximation to 
P(N^n), it might not be reasonable to model the detection process 
in terms of single "look" probabilities. If this is the case, 

Y (t) must be defined by some other means. The visual detection 
model developed in the OEG Report No. 56 is an example of this. 

Generally, in search operations, the target will be moving 
on a relative track C. In some cases, it may be useful to con- 
sider the track to consist of segments c i ,C 2 ' ' * * 'C • Let ^i 
represent the event the target is first detected on track segment 
. The probability of first detecting a target while it moves 
along C can then be expressed as follows: 

n 

P ( r ) = l - exp{- l F(C. ) } 

i=l 

where F(C^) is called the sighting potential on the track seg- 
ment and exp{-F(C^)} = P (T . | H . . . 0 f ^) . It is conven- 

ient to use the concept of sighting potential in the analysis of 
some search types, for example, "ladder searches". 

The track of a target which moves in a plane can be described 
by equations x = x(t) and y = y(t). These equations describe 
the relative track of the target with respect to the searcher 
(often called the observer) if x and y refer to a coordinate 
system in which the searcher is located at the origin. In this 
case, r the range of the target is given by r = [x 2 +y 2 ] >5 . 
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If in a model the detection capability of a searcher 
against a target is assumed to be independent of the time at which 
the target occupies a particular position on the target's track. 



then time can be replaced by target distance along the track, 

ft 

This can be done by using the transformation equation s = 



w(t) dx 



where w is the target's speed relative to the searcher and the 
zero of time has been chosen to coincide with the start of the 
search. If, in addition, the detection capability of the searcher 
against the target is assumed to depend only on the range of the 
target, then the sighting potential for a target's track C for 



a continuous model can be expressed in terms of either t or s 

ft s 

Y [r (t) ] dx = 

0 J 0 



as follows: F(C) = 

dt/ds = l/w(s) . 



Y [r (cr) ] do/w (o) since 
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VII. Experimental Validation of Detection Models 

Consider a detection system which could be satisfactorily 
modeled by a continuous detection rate function that depends only 
on r, the target's range, if a suitable detection rate function 
could be found. Suppose on some basis that a particular detec- 
tion rate function is proposed? An estimate of the suitability 
of the proposed function could be made using standard statistical 
tests if appropriate operational data could be obtained. For 
example, the number of detections at various ranges for straight 
line encounters. If such data were available, its agreement with 
the cumulative distribution function of the target's range at 
detection deduced from the proposed detection rate function could 
be determined by use of the Kolmogorov-Smirnov test. This technique 
could be applied even though the number of encounters in which a 
target is not detected was unknown. In this case, the cumulative 
distribution function to be compared is F R (r|D) where D = 

{detection} and R represents the target's range at detection. 

To derive F R (r|D) or equivalently the density function 
f R (r|D) given a continuous detection rate function, one can pro- 
ceed as follows: First find the joint density function which de- 

scribes the distribution of the target's rectangular coordinates 
(X, Y) at its detection point. To do this, set f x Y (x,y) |D)AxAy = 
f v v (x,y)AxAy/p where p = P(D). This can be done, since 

A , I 

[ (Y=y) flD] = (Y=y) . That is, Y can be considered to take on a 
value only when a detection occurs. Note, p is the value of the 
integral of f v v (x,y) over all pairs of values (x,y) . Now note that 

X , I 
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f (x,y) = . f | v (y | x) f v (x) and that the conditional density 
X / Y Y | A X 

f Y|x(y i x ) can wr; '- tten as 



£ Y | K (y|x) 



( 1/w) Y [r (y ) ] exp{- 



y 



Y [r ( rj ) ] dri/w} . 



— oo 

Given f v (x) and Y(r), f D (r|D) can be found by first trans- 
forming to r and o by using the transformation equations x = r 
cos c and y = r sin a and then finding the marginal distribution 
of R as follows: 



f R (r ) 



2ir 

f r £ (r ,a) da . 

0 



An example of the above procedure will now be given. Suppose 



£ x ( X) = 



l/2a | x | £ a 



x > a 



kh/r 



and y(r) = 



x £ a 



x > a 



The transformation to r and a gives 



r sin a 



(l/2awp) (kh/r 2 )exp{- 



£ R,E (r '°l D) = 



kh(r 2 cos 2 +n 2 ) 



| r cos o | > a 



dn/w} 
r cos a I £ a 



For large a, this can be simplified by letting a °°, after 

f a 

noting that p = l/2a p(x)dx = W/2a where W is the sweep 
width. In the limit, a f R j.(r,o|D) = (kh/wWr 2 ) exp {-kh/wr^ 1-sin a)}, 
and after performing the integration over a, 

f R (r|D) = (W/2r 2 ) {l-$ (kh/wr 2 ) . 

Note, = pf R ( r | D ) which emphasizes the fact that in 

general detection need not take place during a straight line encounter. 
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VIII. Sweep Width Determination for a Random Track Angle 

Now consider the problem of determining the probability p of 
detecting a target during an encounter which could be described 
by a straight line encounter model if the angle (p between the 
target's track and searcher's track were known. Suppose it can 
be assumed that the lateral range and the relative speed 
w = [u 2 +v-2uv cos cp] 2 are values of random variables whose joint 
distribution is known. The probability P of detecting a target 
during a straight line encounter can then be expressed as 



[CO *00 

P = I I p(x,w) f x ^ \|/(X/W) dxdw 
-00 0 

where p(x,w) = p (det | X = x H Y = w) and f x ^(x,w) is the joint 
distribution of the lateral range X and the relative speed Y. 

If X and Y are independent random variables and if 
p(x,w) as ^ function of w is adequately described by a linear 
approximation, then p can be approximated as follows: 



■r 



p (x , w) f x (x)dx 



where 



w = E (<{j ) = 



wfy(w)dw = 



2 ir 



(u 2 +v 2 -2uvcos ip) 2 f^(<p)d<p. 



An average sweep width can be defined by 



W 



pOO 

p(x,w) dx. 

— 00 



If $ is uniformly distributed so f^(<p) = jijp , w is the value 
of an elliptic integral of the second kind. An average sweep width 
based on this value of w should be used in the Random Search 
Formula when either u or v are not appropriate. 
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IX. A Parallel Sweep Search Model 

A model of a parallel sweep or "ladder" search will be 
considered next. In the model, the target's velocity is zero 
and its position is uniformly distributed in a rectangular region. 
The searcher moves along a series of n parallel tracks so that 
n straight line encounters will occur in a complete search of 
the region. The track spacing is s and the tracks are numbered 
from 1 through n. The search geometry is shown in Figure 2. 



b 



1 2 



0 



i 



n 

r 

ns 



* z 



Figure 2. Search geometry for a parallel sweep model. 



If the target's horizontal coordinate is a random variable 

Z in the coordinate system shown in the figure, then P(det|z=z) 
n 

= l-exp[- £ F(x.)] where x. = z - (j-%)s is the target's 

i=l - 1 - 1 

J th 

lateral range from the j — track and 



P (dfet) 



[•ns 

P (det | Z=z)dz/ns. 

J 0 
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In some cases, it is convenient to use the following coordinate 
system: Let I be the random variable which is the number of 

the track closest to the target. Relabel the tracks so that the 

number of the track closest to the target is 0, then = -js + x 

where now j = -K, . . . , 0 , . . . ,M where n = M + K + 1 and the value 

of K and M are determined by the value of I, that is, by 
the location of the target. Now 



m 



P (det | X = x H I=i) = l-exp[- £ F(x-)] = p(x,i) 

j=-k D 



and, since -s/2 < x < s/2, 



P (det 1 1 = i) = 



s/2 



-s/2 



p (x,i)dx/s = p (i) . 



Hence, 



P (det) = l p ( i) P ( I=i) = (1/n) £ p(i) 

i=l i=l 



or 



n 



P (det) = (1/ns) l 

i=l 



s/2 



p (x, i) dx 



-s/2 



If n is large and the detection law is such that y(r) is 

essentially zero for targets at ranges equal to or greater than 

one or two track spacings away, then for almost all values of i 

and j, p(i) will be approximately equal to p ( j ) and it can 

fS/2 

be assumed that p(det) = (n/ns) p(x;s)dx where 

-s/2 

00 

p(x;s) = 1 - exp [- l F (- js+x) ] . 

j=-00 
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As an example, suppose y (r) = kh/r 3 and b is effectively 
infinite. Then F(x^) = 2kh/wx^ . If it is assumed that the 
number of tracks and the track spacing is such that p ( i ) = p ( j ) 
for almost all values of i and j , then 
m 00 

l F (x . ) - £ F(-js+x) = (2kh/w) (tt 2 /s 2 ) esc 2 (ttx/s) 

j=-k 3 j=-°° 

SO 

p(x;s) = 1 - exp{- (2khir 2 /ws 2 ) esc 2 (irx/s) } . 

and 

P(det) = 2${ (tt/ 2) ^ (W/s)} - 1 

h 

where W = 2(2irkh/w) is the searcher's sweep width against the 

target. The ratio W/s is called the coverage factor. A plot of 
P(det) is given as a function of W/s in the National Search 
and Rescue Manual. 

The track length which is expended in an unsuccessful search 
is nb by the assumptions of the parallel sweep model. Since 
the area searched is nsb, the probability of detecting a target 
with a random search in the area with track length nb is given 
by P(det) = 1 - exp{-W/s}. This value could be compared to that 
given by the parallel sweep model in order to obtain a relative 
measure of effectiveness for a parallel sweep search. 
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X . The Optimal Allocation of Search Effort Problem 

The problem of determining an optimal allocation of effort 
for a search is discussed next. Here, an allocation of effort is 
an optimal allocation if it maximizes the probability of detecting 
the target. 

The class of searches which will be considered are those 

which can be described by a random search model that satisfies the 

following conditions: The target is fixed in a region of area A. 

The region consists of m subregions which are determined by the 

th 

conditions that in the i— subregion the sweep width is a con- 
stant and the probability density function is a constant 

= p^/A^ where A^ is the area of the subregion and p^ is 
the target's prior probability of being there. (The subregions 

are assumed to be numbered so that W, P. > W~P~ > . .. > W P .) If a 

1 ± 2. 2. n n 

target is detected in a subregion, the searcher must be in the 
subregion when it is detected. The probability of detecting the 
target is then 

n 

P(det) = \ {l-exp(-(pj ) }p- A. 

i=l 

where cp^ = W^£^/A^ is called the search effort density. 

Now consider the problem of determining the solution set 

cp . which maximizes P(det) subject to the constraints 

1 °P n 

SL = £ (A . /W . ) <p . and <p . £ 0. Note, A is the available track 

i=l 

length. This is a nonlinear programming problem for which the 
solution set is: 
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In W^p i - (1/A) £ (Aj/Wj) In Wj Pj + Jt/X i in ft 



ft 



<P 



1 op 



not in ft 



where X = £ A./W. and ft = {l,2,...,k} where k is determined 
ft ^ ^ 

by the condition that if k + 1 were included in ft then <p^ + ^ < 
Suppose the model described above were used to determine 
an optimal allocation of track length for a search. The result 
would be a set of values l, 



'i op 



= (A^/VT) cp^ o p which maximize 



the probability of detecting the target. Given detection will 
occur, in what order should the subregions be searched in order 
to minimize the expected time until the target is detected? The 
following procedure will effectively minimize the expected track 
length to detection and gives a suitable order assuming the 
searcher's speed remains constant so that minimizing the expected 
track length to detection is equivalent to minimizing the expected 
time to detection. 

Divide H into minimum increments AJl consistent with 
the random search model and then allocate the increments in the 
following order: Allocate A SL to the Region 1. This will be 

assumed to be effectively consistent with the prescription deter- 
mined by the formula for <p^ given AJl. If this search is 

unsuccessful, determine the optimal allocation prescribed for 
2A£. This will effectively indicate whether the next increment 
AJl should be assigned to Region 1 or Region 2. Continue in this 
fashion until detection or until Jl has been expended. Note, 
this is an approximate procedure in the sense of the model. 



26 



For a discussion of the above procedure, let E(l|l^£) 
symbolize the expected track length to detection given detection 
will occur by the time track length £ has been used. That the 
above search procedure will effectively minimize E(l|l££) for 
an optimal allocation of track length is suggested by the follow- 
ing argument: Divide £ into segments A£ as indicated above, 
then the probability that the target is detected on or before the 
i— step of the search P(L^iA£) is, in effect, a maximum, 
since no other allocation of search effort iA£ will give a 
significantly greater value for any value of i. Hence, 

F (iA£|L^£) = P(L^iA£|L^£) is also effectively a maximum for any 
value of i since P(L^iA£|L^£) = P (L^iA£) /p (L^£) and P(L££) 

must be equal to its maximum value for any alternate procedure. 

n 

Therefore, E(l|l^£) = \ [1-F T (iA£ | L^£) ] is effectively a 

i=l L 

minimum for the procedure. 

The above procedure is the one which would have been followed 
if after each unsuccessful search one calculated new prior proba- 
bilities using Bayes procedure and then did the search prescribed 
for a track length A£ with the new priors. For a discussion of 
this point, see OEG Report No. 56. An application, in effect, is 
discussed in Reference 4. 
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XI. Target Position Probability Distributions Which Change in Time 
Consider the following model in which a target moves with 
constant course and speed in a plane with its position specified 
in a rectangular coordinate system: The target's course and speed 

are independent random variables $ and U respectively 
whose density functions and f^(u) are 9i- ven * The joint 

density function f v v (x,y;t) of the target's coordinates is 

A , I 

needed, but only f v v (x,y;0) the joint density function for time 
zero is given. 

For such cases, f v v (x,y;t) can be found in principle as 
follows: First it is convenient to use polar coordinates and with 

x = p sin a and y = p cos a to set 

f ( p , a ; t ) = f X j(P sina,pcosa;t). 

Note, to first order, P{target is in AA at t} = f(p,a;t)AA. 

North 




Figure 3 . Problem geometry . 
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Referring to Figure 3 



f<*> f 2 TT 



f (p,a;t) AA = 



0 0 



f [r (u, <}>) , 9 (u,<j>) ; 0] AA' f ^ (<f>) f y (u) d<J>du 



where the functions r(u,<}>) and 0(u,<f>) are such that if the 
target's polar coordinates were r(u,<J>) and 0(u,<j>) at time 0, 
then with the speed and course u and <j> they would be p and a 
at t. They are determined by r cos 6 = p cos a - ut cos and 
r sin 0 = p sin a — ut sin <j>f so r 2 = p 2 + (ut) 2 - 2p ut cos ( a— <j) ) . 
The integration is, in effect, a sum over all possible pairs of 
values of u and <j> . Each pair determines a vector ut which 
translates target positions at time 0 to target positions at 
time t such that Aa is AA' translated without rotation or 
distortion. So, AA' = AA, independent of u and <J> . Therefore, 



f ( p , a ; t ) = J j f [r (u,()>) ,0 (u,<{>) ; 0] f $ ((|>) fytujd^du 



0 0 



As an example, suppose 





2a 



with f $(<}>) = / 0 < <f> £ 2ir, but with U = u known. In this 



case , 



f [r (u,<j>) ,0 (u,<f>) ; 0] = 2vo r e 




and, using r 2 = p 2 + (ut) 2 - 2p ut cos(a-<J>) 



29 



f (p,a; t) 



2ira 2 j 



2 it - -~r [ p 2 + ( ut ) 2 — 2 p ut cos(a-t)) ) ] 



d({) 

2tt 



3 _ - [p 2 + (ut) 2 ] 1 ,2 tt ^ cos(a-<f>) 



2ttc 



t e 



— \ 
2 IT J 



d<t> 



X rr [p 2 +(ut) 2 ] 



2ira 



■r e 



2a" tr " ' """ J T riPut'i 
J 0 



where J 0 is the zero order Bessel function of the first kind 
and i = /-l. To obtain f v (x,y;t), substitute P 2 = x 2 + y 2 

A / x 

in the above result. In this case, since f x , Y (x,y;t) is symmet- 
ric about the origin, the marginal distribution of the target's 
range R from the origin is of interest. In general, 



f H Q (P,a;t) = Pf x y (p sin a,p cos a ; t) = pf(p,a;t) 



and 



r 2 it 



f R (p; t) = 



f R,e (pf 



so, in this case. 



f R (P;t)=^e 20 



1 T [p 2 +(Ut) 2 ] 






The distribution function F R (p;t) for this case is given in 
Reference 5. 
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